LONG TIME, LARGE SCALE LIMIT OF THE WIGNER 
TRANSFORM FOR A SYSTEM OF LINEAR 
OSCILLATORS IN ONE DIMENSION 
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Abstract. Wc consider the long time, large scale behavior of the 
Wigner transform W e (t, x, k) of the wave function corresponding to 
a discrete wave equation on a 1-d integer lattice, with a weak multi- 
plicative noise. This model has been introduced in [2] to describe a 
system of interacting linear oscillators with a weak noise that con- 
serves locally the kinetic energy and the momentum. The kinetic 
limit for the Wigner transform has been shown in [J] . In the present 
paper we prove that in the unpinned case there exists 70 > such 
that for any 7 € (0,70] the weak limit of W e (t/e 3 / 2 t , x/e 7 , k), 
as e -C 1, satisfies a one dimensional fractional heat equation 
d t W{t,x) = -ci-dlf^Wit^x) with c > 0. In the pinned case 
an analogous result can be claimed for W e (t/e 27 , x/e 1 , k) but the 
limit satisfies then the usual heat equation. 



1. Introduction 

In the present paper we are concerned with the asymptotic behavior 
of the Wigner transform of the wave function corresponding to a dis- 
crete wave equation on a one dimensional integer lattice with a weak 
multiplicative noise, see (12. ip below. This kind of an equation arises 
naturally while considering a stochastically perturbed chain of oscil- 
lators with harmonic interactions, see [2] and also [IE]. It has been 
argued in [2] that, due to the presence of the noise conserving both the 
energy and the momentum (in fact the latter property is crucial), in 
the low dimensions (d = 1, or 2) conductivity of this explicitly solvable 
model diverges as iV 1 / 2 in dimension d — 1 and logiV, when d — 2, 
where N is the length of the chain. This complies with numerical re- 
sults concerning some anharmonic chains with no noise, see e.g. [17] 
and [16] . We refer an interested reader to the review papers [7J [16] and 
the references therein for more background information on the subject 
of heat transport in anharmonic crystals. 

It has been shown in [I] that in the weakly coupled case, i.e. when the 
coupling parameter e is small, the asymptotics of the Wigner function 
W e (t,x,k) (defined below by (12. 8p . with 7 = 0), that describes the 
resolution of the energy in spatial and momentum coordinates (x, k) 
at time t ~ e _1 , is given by a linear Boltzmann equation, see (12. 9 p 
below. Furthermore, since in the dimension d — 1 the scattering rate 
of a phonon is of order k 2 for small wavenumber k, in the unpinned 
case (then the dispersion relation satisfies u'(k) signfc, for \k\ <C 1) 

the long time, large space asymptotics of the solution of the transport 
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equation can be described by a fractional (in space) heat equation 

d t W(t,x) = -c(-d 2 x ) 3/A W(t,x) 

for some c > 0. The initial condition W(0,x) is the limit of the average 
of the initial Wigner transform over the wavenumbers, see [13] and also 
[31 [TJ5]. Note that the above equation is invariant under time-space 
scalings t ~ t'/e 37 / 2 , x ~ x'/e 7 for an arbitrary 7 > 0. This suggests 
that the fractional heat equation is the limit of the Wigner transform 
in the above time space scaling. In our first main result, see part 1) of 
Theorem 12. II below, we prove that it is indeed the case when 7 £ (0, 70] 
for some 70 > 0. 

On the other hand, in the pinned case, i.e. when the dispersion 
relation satisfies u'(k) ps 0, as \k\ <C 1, one can show that the solution 
of the Boltzmann equation approximates the regular heat equation 

d t W(t,x) =cd 2 x W{t,x). (1.1) 

In part 2) of Theorem 12.11 below we assert that if the Wigner transform 
is considered under the scaling (£, x) ~ (f '/e 27 ', x'/e 1 ), for 7 £ (0, 70] and 
some 70 > 0, then it converges to W(t, x), as e <§C 1. The coefficient c in 
the heat equation ( II. ip differs from the thermal conductivity coefficient 
calculated explicitly in the pinned case with the help of the Green-Kubo 
formula in [TJ, see Theorem 1. In fact, its computation, see formulas 
( I6.12p - ( I6.13P below, requires solving a Poisson equation (I6.14p and 
does not lead to an explicit formula. Finally, we mention also the results 
concerning the diffusive limits for the Wigner transform of a solution 
of the wave equation on a lattice with a random local velocity in the 
weak coupling regime (see [ISj), for the geometric optics regime for the 
wave equation in continuum (see [15]) and in the case of Schrodinger 
equation in the radiative transport regime (see [12]). 

The proof of Theorem 12.11 is made of two principal ingredients: the 
estimates of the convergence rate for the Wigner transform towards the 
solution of the kinetic equation, see Theorem 13.2 1 below in the unpinned 
case (resp. Theorem 13.51 for the pinned case), and the respective rate of 
convergence estimates for the solutions of the scaled kinetic equation, 
see Theorem 13.31 (resp. Theorem 13. 6p . To prove the latter we show 
two probabilistic results: Theorems 15.51 and 15.81 that are of interest on 
their own. They provide estimates of the rate of convergence of the 
characteristic functions corresponding to a scaled additive functional 
of a stationary Markov chain towards the characteristic function of an 
appropriate stable limit and the respective result in the continuous time 
case. 
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2. Description of the model and preliminaries 

2.1. Discrete wave equation with a noise. We consider a discrete 
wave equation with the multiplicative noise on a one dimensional inte- 
ger lattice, see [2], 



(2.1) 



dpx_ _ _ g , 
^ dt 



Here (p,q) = {(p x ,q x ), x G z}, where the component labelled by x 
corresponds to the one dimensional momentum p x and position q x . The 
Hamiltonian corresponds to an infinite chain of harmonic oscillators 
and is given by 

y ez y,y'ez 

The interaction potential {a x , x G z} will be further specified later 
on. The noises {£x(t), x G z} are defined by the following stochastic 
differentials 

d£\t) = y/i Yl (Y x+k p x )odw x+k {t), (2.2) 

fc=-l,0,l 

understood in the Stratonovich sense. Here 

Y x := (p x - px+iidp^ + (p x +i - px-i)d Vx + (p x _i - p x )d Px+1 

and {w x (t), t > 0}, x G Z are i.i.d. standard, one dimensional Brow- 
nian motions over a certain probability space (Q, J 7 , P). Note that the 
vector field y x is tangent to the surfaces 

p*_i + P* + Px+i = const (2.3) 

and 

p^-i + p x + p x+1 = const (2.4) 

therefore the system (12. ip conserves the total energy and momentum. 
System (12. ip can be rewritten formally in the Ito form: 

dq y (t) = p y (t)dt (2.5) 

dpyit) = [-(a * q(t)) y - * p (*))„] dt, 
+v / e XI ? iH-fcP»(*)) d ' u, iH-fc(*)> 2/ ez - 
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Here /3 y := A/3£ 0) , with 



Pf ] = { -l, y = ±1 



4, y = 

1, = ± 
0, if otherwise. 



The lattice Laplacian for a given g : Z — >■ C is defined as Ag y := 
9y+i + fl'jz-i — 2# r 

2.2. Formulation of the main results. To describe the distribution 
of the energy of the chain over the position and momentum coordi- 
nates it is convenient to consider the Wigner transform of the wave 
function corresponding to the chain. Adjusting the time variable to 
the macroscopic scale it is defined as 

^)(t) :=a)*q +ip Qj. (2.6) 

Here u is the inverse Fourier transform, see (12.171) . of the dispersion 
relation function given by ou(k) = \/ 6t(k), with a[k) the direct Fourier 
transform of the potential, defined on T- the one dimensional torus, see 
(12.16p . Suppose that the initial condition in (12.1 ft is random, indepen- 
dent of the realizations of the noise and such that for some 7 > 

limsupV e 1+ m^ v e \0)\ 2 } e < +00. (2.7) 

Here (-) e denotes the average with respect to the probability measure 
fi e corresponding to the randomness in the initial data. In fact, since 
the total energy of the system ^2 y€Z |V4^WI 2 i s conserved in time, see 
Section 2 of [4], an analogue of condition ( 12. 7p holds for any t > 0. 

The (averaged) Wigner transform of the wave function, see pE], is a 
distribution defined as follows 

(W £ ^t), J) := — I e i2 <*'-^ k J* ^— (x + *'), kj 

^\t)Y^(t)] (2.8) 



for any J belonging to S - the Schwartz class of functions on R x T, see 
Section [231 Here E e is the average with respect to the product measure 
/4<g)P. It has been shown in [I], see Theorem 5, that, under appropriate 
assumptions on the potential «(•), see conditions al) and a2) below, the 
respective Wigner transforms W e ${t) converge in a distribution sense, 
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as e — y 0+, to the solution of the linear kinetic equation 

d t U(t, x, k) + -±±d x U(t, x, k) = CU(t, x, k), (2.9) 
2ir 

where £ is the scattering operator defined in ( I2.35p . Our principal re- 
sult concerning this model deals with the limit of the Wigner transform 
in the longer time scales, i.e. when 7 > 0. It is a direct consequence 
of Theorems 13.11 and 13.41 formulated below that contain also the infor- 
mation on the convergence rates. Before its statement let us recall the 
notion of a solution of a fractional heat equation. Assume that Wo is 
a function from the Schwartz class on R. The solution of the Cauchy 
problem for the fractional heat equation 

dtW&x) = -J^Y b (-9l) b/2 W(t,x) (2.10) 

with b, c > and W(0, x) = Wo(x) is given by 

W(t,x) := f e i27rxp -W bt W (p) dp, 
Jr 

where 

W (p) := [ e- i2wxp W (x) dx 
Jr 

is the Fourier transform of Wq(x). 

Theorem 2.1. Suppose that potential {a y , y G z} satisfy assumptions 
al) — a2) formulated in Section [Ol Then, the following are true. 

1) Assume that d(0) = (no pinning case), 7 G (0,2a/3) for 
some a G (0, 1] and 

K aa := hmsupe 1 ^ / (|#(fc)| 2 ) e ^ < +00, (2.11) 

where ^^{k) is the Fourier transform of the initial condition 
i/jx\o). Suppose also that for some W with the norm (I2.24p 
finite we have 

lim (W e7 (0), J) = [ W (x,k)J*(x,k)dxdk, (2.12) 

' J RXT 



for all J G S. Then, for any t > we have 



hm (w ea \-£fiJ , J) = J W(t,x)J*(x,k)dxdk, (2.13) 
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where W(t,x) satisfies (12.101) with b = 3/2 and the initial 
condition given by 

W{0, x) := W {x, k)dk. (2.14) 

The coefficient c is given by (13. 5p . 
2) Suppose that a(0) > (pinned case), 7 G (0,1/2) and condi- 
tions (I2.1ip and fl2. 12j) for Wq with the norm (I2.24p finite are 
satisfied. Then, for any J G S and t > we have 

Urn /w e>7 (J^J , Pj = J W(t,x)J*(x,k)dxdk,, (2.15) 

where W(t,x) is the solution of the ordinary heat equation, i.e. 
(I2.10p with b = 2, with the initial condition given by (I2.14p and 
coefficient c as in (16 . 1 2[) . 

2.3. Fourier transform of the wave function. The one dimensional 
torus T is understood here as the interval [—1/2, 1/2] with its endpoints 
identified. Let e r (k) := exp{— i2irrk}, r G Z. It is a standard orthonor- 
mal base in L^.(l) - the space of complex valued, square integrable 
functions. The Fourier transform of a given square integrable sequence 
of complex numbers {g y , y G z} is defined as 

g{k) = Y J 9 v e y {k), keT (2.16) 

and the inverse transform is given by 

fv = Jel{k)f{k)dk, yez (2.17) 

for any / belonging to L^(t). 

A straightforward calculation shows that ip^ (t, k) - the Fourier trans- 
form of the (complex valued) wave function given by (12.61) - satisfies 
the following Ito stochastic differential equation, cf. formula (7.0.7) of 

i, 

d^ e \t) = A e [^\t)]dt + ^Q[^\t)}(e r )dw r (t), (2.18) 
^)(0) = ^ , 

where ^0 is the Fourier transform of ^^(0). Here, a (nonlinear) map- 
ping A : L\ (T) -> L\ (T) is given by 

A e [f)(k) := —-u(k)f(k) - ffi a ^ E L '(T), (2.19) 
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where 



/+(*):=/(*) and /_(*):=/*("*), 



£(fc) = 8sin 2 (vrfc) [l + 2 cos 2 (tt£;)] . 



(2.20) 



For any g £ L 2 (l) a linear mapping : L 2 (t) — > L 2 (t) is given by 



r(Jfe, fc') := sin(27rA;) + sin[2vr(£; - fc')] + sin[27r(A;' - 2Jfe)] 
= 4 sin(7r/c) sin[7r(fc - k')] sin [(2fc - fc')7r] , k, k' £ T. 

Finally, {w r (t), t > 0}, r £ Z are i.i.d. one dimensional, standard 
Brownian motions, on a probability space (Jl,J-,P), that are non- 
anticipative w.r.t. the given filtration {J-" f , t > 0}. 
We assume, as in [1], that 

al) {a y , y £ z} is real valued and there exists C > such that 

K| < Ce-l y l/° for all y £ Z, 
a2) a(/c) is also real valued, a(k) > for k ^ and in case 

a(0) = we have d"(0) > 0. 

The above assumptions imply that both y H- a y and k H- d(/c) are real 
valued, even functions. In addition a £ C°°(t) and if a(0) = then 
a(k) = k 2 <f)(k 2 ) for some strictly positive £ C°°(t). This in particular 
implies that, in the latter case, the dispersion relation u(k) = \J 'a(k) 
belongs to l7°°(t \ {0}). 

It can be easily checked that under the hypotheses made about the 
potential a y the mapping given by fl2~T9l is Lipschitz from L 2 c (t) to 
itself and £ rez \\Q[g}(e r )\\ 2 L2 < C\\ g\\ 2 L2 for some C > and all g £ 
L 2 (t) so Q[g] is Hilbert-Schmidt. Using Theorem 7.4, p. 186, of [9] one 
can show that for any L 2 (T)-valued, J-o _measura ble, initial data 
there exists a unique solution to (I2.18p understood as L 2 (t) - valued, 
continuous trajectory, adapted process {i/j( e '(t), t > 0} a.s. satisfying 
(12.18)) . In addition, see Section 2 of [1], for every initial data ip^ £ 
L 2 (t) we have 




where 



||^( e )(t)|| L 2 = const, Vt>0, P a.s. 



(2.21) 



2.4. Wigner transform. 
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Some function spaces. Denote by S the set of functions J : R x T — > C 
that are of C°° class and such that for any integers I, m, n we have 
su Pp,fc(l + P 2 ) n ^ 2 \9 l p d'^ L J(p, k)\ < +00. For any a G R we introduce the 
norm 

||J|k := / (l+p 2 ) a/2 sn V \J(p : k)\dp : (2.22) 

JR k 

for a given J G S. By A' a we denote the completion of S in the norm 
|| ■ H^. Note that A' a is dual to A a defined as the completion of S in 
the norm 

||J|U :=sup(l+p 2 )- a / 2 / \J(p,k)\dk. (2.23) 
p Jt 

We use a shorthand notation A := Ao and A' := A' . With some 
abuse of notation by (■, •) we denote the scalar product in L^,(t) and 
the extension of 

(Ji,J 2 ):=/ Ji(p,k)Jl(p,k)dpdk 
Jrxt 

from S x S to A x A'. 

We shall also use the space £y& obtained by completion of S in the 
norm 

||J||^:=sup(l+ P r 2 /^#^ (2-24) 
p Jt |"m 

When 6 = we shall write B a instead of B a fl. 

Random and average Wigner transform. For a given e > let i/j^(t) be 
a solution of (12.181) with a random initial condition ijj^(0) distributed 
according to a probability measure /i e on L 2 (t). Define 

W e (t, p, k) := ( (>)) * (t,k - |) ^ « (f , fc + |) ^ (2.25) 

and 

? e (f , p, fc) := (t, -k + |) # e > (t, * + f ) \ , (2-26) 

where, as we recall, (-)<= is the average with respect to the initial con- 
dition. Using (P22D and (12T2TT) we conclude that both W e (t) and Y e (t) 
belong to L 1 (P; A) - the space of .4- valued random elements possessing 
the absolute moment, for any t > 0. We also introduce the average ob- 
jects W t (t,p, k) and Y e (t,p, k) using formulas analogous to (I2.25P and 
(I2.26P with (-) e replaced by E e corresponding to the average over both 
the initial data and realization of Brownian motion. 
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The (averaged) Wigner transform W ea {t) is defined as 



(W^(t),J):-- 

where J G A' and 

J(x, k) : 



r l+7 



W e (t,e' ( p,k)J*(p,k)dpdk (2.27) 



exp {i2irpx} J(p, k)dp. 



The anti-transform Y t ^{t) is defined by an analogous formula, with W e 
replaced by Y e . 

2.5. Evolution of the Wigner transform. Using Ito formula for 
the solution of (12.181) . see Theorem 4.17 of [9], we conclude that 

(2.28) 



dW e (t,p, k) 



+ E ( (*> k - f ) (t, k + 1 ) } ) dt 

+dM { t e \p,k), 

where the {-M| e) , t > 0} is an {J^t, t > 0}-adapted local martingale, 
given by 

Mf{jp,k) (2.29) 
:= E T ((« (e) (*)* (* - |) ^ (s, k + | )} d Wj (s) 

+ e r («> * - 1 ) (* + y ) ) e 

In order to guarantee that the stochastic integrals defined above are 
martingales, not merely local ones, we need to make an additional 
assumption that fi t has the 4-th absolute moment. Taking the expec- 
tation of both sides of (12.281) with respect to the realizations of the 
Brownian motion we conclude that W e (t,p, k) = EW e (t,p, k) satisfies 



d t W e (t,p, k) 



i5 e u(p, k) + -p e (p, k) 



4 / p e (k,k- k',p)E (p {e) T U,k' -^-)p^ (t,k'+^-) dk', (2.30) 



w e (t, P ,k) + nf\t,p,k) 



cp 
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P [e \t,k) :=-[^\t,k)-(^y(t,-k) 

( ^P \ ( ^P 
Pe(k, k',p) :=r[k- —, k'J r[k + —, k' 

6 £ u(p, k) : 



€ 

1 

2 L 



oj [k + 
$ (k + 



cp 
~2 
ep 

y 



oj k 



ep\ 
2 J 



+ (3{k- 



ep 



(2.31) 



and 



H?\t,p, k) := i {$ (k - |) F e (t,p, fc) + ^ (k + |) <(*, -p, k) 

with Y t (t,p, k) = EY t (t,p, k). 

Formula ( I2.30p remains valid also when only the second absolute mo- 
ment exists. This can be easily argued by an approximation of the ini- 
tial condition by random elements that are deterministically bounded. 

Since the momentum, that is the inverse Fourier transform of p^ e \t, k), 
is real valued, the expression under the expectation appearing on the 
right hand side of (I2.30p is an even function of k', thus the last term 
appearing on the right hand side of the equation can be replaced by 



4 / R e (p,k,k')E e (pU)* U,k' -^)pV (t,k'+^-) dk', (2.32) 



cp 



where 



Note that 



R e (p, k, k') := - ^2 Pe{K k + ik',p). 



t=±i 



R(k, k') := R (p, k, k') = - [r 2 (k, k - k') + r 2 (k, k + k' 

= 8 sin 2 (7rfc) sin 2 (vrA;') {sin 2 [n(k + k')} + sin 2 [n(k - k')]} . 
The following relation holds 

4 / R(k,k')dk' = j3(k), VfcGT. (2.33) 



We conclude therefore that W e (t,p, k) satisfies the following 

(d t W e (t),J) = (W e (t),(iB + £)J) + (1l e (t),J), VJeS. (2.34) 

where Bf(p,k) := poj'(k)f(p,k), for any f E S. We let u/(0) := in 
case oj is not differentiable at 0. In addition, 

C : = £ (0) , (2.35) 
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where for each e G [0, 1] operator acts on S according to the formula 
C^f(p, k):=2^ R £ (p, k, k')f(p, k')dk' - \Hk)f{p, k) 

= 2 ^R e (p,k,k')[f(p,k')- f(p,k)]dk', feS, (2.36) 

and extends to a bounded operator on either A, or A'. Finally, 

K e (t,p,k) :=nU(t,p,k)+K?(t,p,k), (2.37) 

with 

TZ^(t,p, k) := [i (pco'(k) - 6 e u(p, k)) + (jrM _ £)] w e (t,p, k), 



nf\t, P ,k) :=nf\t, P ,k)+ Tz^(t,-p,k) 



(2)/ 



(2.38) 



where 

K?\t,p,k) := 1/3(A;-| W 



Y e (t,p,k)- / R e (p,k,k')Y e (t,p,k')dk'. 



Similar calculations show that 

2i 



jY.(t, P ,k) 



u e (p, k)Y e (t,p,k)+U e (t,p,k), 



(2.39) 



where 



x 



e(t,p,k) := --P e (p,k)Y e (t,p,k) - <Tl<T2 [ Pe(p,k,k') 

(T1,<7 2 =±1 ^ T 

Ee {^) (f , fc - A/ - I ) (</©* (t, fc + fc' + I) } dk' 



and 



w ( k + — ) + u) [ k 



uj t (p,k) :-- 



p e (p, k, k!) :— r [k — — , k \ r [k 



ep 
ep 



—k' 



2.6. Probabilistic interpretation of the kinetic linear equation. 

Denote by U(t,p, k) the Fourier transform of the solution of (12.91) in 
the x variable. Let K t (k) be a T- valued, Markov jump process, defined 
over (Q, J 7 , P), starting at k, with the generator C. Suppose also that 
U G A. Then 

d t U(t) - iBU{t) = CU(t), 17(0) = U (2.40) 
is understood as a continuous A- valued function U (t) such that 

(U(t), J) - (U , J) = [ (U(s), (iB + C) J)ds (2.41) 



12 



TOMASZ KOMOROWSKI AND LUKASZ STEPIEN 



for all J £ A'. It is well known that this solution admits the following 
probabilistic representation 



exp -ip / u' (K 8 (k)) ds } U (p, K t (k)) 



U(t,p, k) := E 



Here E is the expectation over P. For given J £ A' we let 



(2.42) 



J(t,p, k) 



E exp u' (K.(k)) dsj J{p,K t {k)) . (2.43) 

Therefore J(t) £ A'. Using the reversibility of the Lebesgue measure 
under the dynamics of Kt, we conclude that the law of {K t - S , s £ [0, £]} 
and that of {K s , s £ [0,t]} coincide. Hence, 

(U(t),J) = (U Q ,J(t)). (2.44) 

Likewise, using the definition ofH £ (i) (see (12.371) ). from (12.341) and the 
Duhamel formula we get 

(W e (t)i J) = (WM J(t)) + f(Ke(s), J{t - s))ds, VJ £ S. (2.45) 



3. Convergence of the Wigner transform 
For a given a £ R define the norm 

\f(k)\ 2 



2a 



-dk 



1/2 



Recall also that /x e is the distribution of the initial data for equation 
(F2~T8|) . We assume that: 

A ) for a given a > 

/C a , 7 := limsupe 1+ ^ / \\f\\ 2 H - a ^(df) < +oo. (3.1) 

Let := A^o,7- 

3.1. No pinning. Since 



J 



r l+7 



RXT 



3 7 /2 ' 



pe 7 , J*(p, k)dpdk, 



part 1) of Theorem 12. II is a consequence of the following result. 

Theorem 3.1. Assume that to > 0, d(0) = and a £ (0,1] is such 
that (13. ip holds. Then, for any 7 £ (0, 2a/3), 

" 3 



< 7' < 7 mm 



1 



13' a+ 1 



(3.2) 
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and b > 1 one can find C > such that 
; 1+7 ^ ( t 



RXT . 



-w f 



.37/2 ' 



pe 7 , fc - W (p)e~ 



e|p| 3 / 2 t 



J*(p, k)dpdk 



A+7 



< 

+C , e 2a ||J| 



■W e (0,pe\k)-W (p, k) 



-c\p\ 3 / 2 t T* 



J* (p) dpdk 



\W< 



0\\B aj 



t 



\J\ 



|J|| B Je 7 '+C7C 7 e||J|U, + 



A' 



2a+l 



K. 



7 .3 7 /2 



1 

(3.3) 



for all e G (0, 1], ji t , t > t , J £ A' 5 H B afi and W G B a . Here 



W (p) 



and 



Wq (p, k) dk and J{p) = 

3/4 



J (p, k) dk 



n 2 a" 



(0)V 



(3.4) 



(3.5) 



Proof. Denote by U e (t,p, k) the solution of (12.41j) with the initial con- 
dition U e (Q,p,k) = Woipe^.k). The left hand side of (I3.3P can be 
estimated by 



r l+7 



r 3 7 /2 



pe 7 , k 



r 3 7 /2 



,pe 7 ,/c 



J*(p, k)dpdk 



,37/2 ■ 



pe 7 , k - W (p)e~ 



i|p| 3 / 2 t 



J*(p, k)dpdk 



(3.6) 



Denote the terms appearing above by J\ and J 2 respectively. The proof 
is made of two principal ingredients: the estimates of the convergence 
rates of the averaged Wigner transform of the wave function given by 
(I2.18P to the solution of the linear equation (I2.40p (that correspond to 
the estimates of J\) and further estimates for the long time-large scale 
asymptotics of these solutions (these will allow us to estimate J-i). To 
deal with the first issue we formulate the following. 

Theorem 3.2. Suppose that {^ (e) (t), t > 0} is the solution of f!2TTgj) 
with coefficients satisfying al)-a2) with d(0) = and A a ) for some 
a G (0,1]. Assume also that U(t) and J(t) are given by (12.421) and 
(I2.43P respectively. Then, there exists a constant C > such that 



A+7 



■W e (t)-U(t),J 



< 



,1+7 



-W e (0)-U o ,J(t) 



+CeJC^t + 1)|| J\\ A[ + C e'HitlC, + /C a , 7 )|| J\\ A 



2a+l ' 



(3.7) 
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for all J eS, ee (0, 1] and t > 0. 



We postpone the proof of this theorem until Section H~Tl proceeding 
instead with estimates of J\. Using (I3.7H we can write that 



Ji< 



dpdk(5W) t (p, k) 



(3. 



XE 



J* (p, K t/e ^ /2 (k)) exp I -ipe 1 J u'(K s (k))ds 



-CJCj€\\ JelU'j 



r 3 7 /2 



Ce 2a \\J elu , 



t 



2a + l f 37/2 \ 7 3 7 /2 



/C- 



,7 / ' 



where J e (p,k) := e 7 J(pe 7 , k) and 



(<W) e (p,A;) := 



r l+7 



-W e (0,pe\k)-W (p, fc). 



Since || J e |U' < II^IU' f° r every e G (0, 1] and b > 0, we obtain the last 
two terms on the right hand side of (13.81) account for the last two terms 
on the right hand side of ( 13. 3ft . 

Denote the first term on the right hand side of (13. 8p by X. We can 
write that X < T\ + X 2 , where 



Xi 
X 2 



(<JW) £ (p, fc)e- 



-c\p\ 3 / 2 t 



J* (p) dpdk 



RXT 



(SW) e (p,k) \e-^ i/2t r (p) 



-E 



-ipe 1 J u'(K s (k))ds 

Here for any function / : R x T — > C we denote 

7(p):= / /(p,fc)dfc. 



dpdk 



Term Xi accounts for the first term on the right hand side of (13. 3p . 
Using the reversibility of the Lebesgue measure under the dynamics of 
K t (see ( I2.44p ). we conclude that 



J* (p,k)\e-W 5/2t (6W) e (p) 



-E 



(SW) e (p, Kt/^/iik)) exp I -ipe 



' u'{K s {k))ds 



dpdk 
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To estimate I2 (and then further to estimate J2) we need a bound on 
the convergence rate of the scaled functionals of the form (I2.42p . Let 

W(t,p) := W (p)exp{-c\p\ 3 / 2 t}. (3.9) 



Theorem 3.3. For any t > 0, a G (0, 1], b > 1 and 7' as in (13. 2 j) 

tnere exists a constant C > swcn inai 



RXT 



jw(t,p)-E 



Wo (p,i^ i/e 3 7 / 2 (/c)) exp 



w'(^(A;))d 



x J*(p, k)dpdk 



<ct\\w \\ B M\J\U + \\Jh a ,y 



for all e G (0, 1], t > W G B a and J G A' 5 D £ a , fe . 

The proof of this result shall be presented in Section 16.11 Using the 
above theorem we can estimate 

X 2 <Ct(||J|U + ||J|kJe^sup 



r l+ 7 



\W e (0,pe\k)\ + \W (p,k) 

p Jt l ^ 

<Ct(/C 7 , a +||Wo|k)(||J|U 5 + ||J|k ; Je^. 

Invoking again Proposition I3.3[ this time to estimate J2, we obtain 
that 

^2<Ct||M/o|| Sa (||J||^ + ||«/|kJ e7 '- 
The above estimates account for the second term on the right hand 
side of (13.31) . thus concluding the proof of the estimate in (13. 3ft . □ 

3.2. Pinned case. Part 2) of Theorem 12.11 is a direct consequence of 
the following result. 

Theorem 3.4. Assume that a(0) > and t > 0. Then, for any 
7 G (0, 1/2), a G (0, 1] and 



< 7' < 



CZ7 



(3.1i; 



one can find C > such that 



RXT 



r l+7 



( t 



-w 



r27 



pe^,k) -W {p)e- dp2t 



J*(p, k)dpdk 



< 



RXT 



r l+7 



■W e (0,pe\k)-W (p, k) 



e -cp 2 tj* ^ dpdk 



+CK,^t\\ J\\ A + Ct{K^ a + ||W |k)(|| AU + \\ J hJ^' 

(3.12) 

for all e G (0, 1], t > t , J G A!± n i3 a , 6 and W e B a . i/ere W (p), J(p) 
and c are respectively given by H3.4[) ana 1 ( 16 . 1 2[) below. 
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Proof. We proceed in the same fashion as in the proof of Theorem 13.11 
so we only outline the main points of the argument. First, we estimate 
the left hand side of (I3.12p by an expression corresponding to (13. 6p . 
The first term is estimated by an analogue of Theorem 13.21 that in this 
case can be formulated as follows. 

Theorem 3.5. Assume that conditions al)-a2) and Aq) hold. In ad- 
dition, we let a(0) > 0. Then, the average Wigner transform W e (t) 
satisfies the following: there exists a constant C > such that 



r l+7 



■W e (t) — U(t),J 



< 



r l+7 



■W e (0)-U o ,J(t) 



Ce\\ Jll^/Cyt 
(3.13) 



for all e E (0, 1] and t>0. 

The proof of this result is presented in Section 14.21 In the next step 
we can estimate the rate of convergence of the functional appearing in 
the formula for the probabilistic solution of the linear Boltzmann equa- 
tion towards the solution of the heat equation given in the following 
theorem. Let 

W(t,p) :=W (p)exp{-cp 2 t}. (3.14) 

Theorem 3.6. For any t > 0, a E (0, 1], b > 1 and 7' as in (I3.1ip 

there exists a constant C > such that 

W (p, K t/e2j {k)) exp <( -%pe< / u\K 8 {k))ds 



\W(t,p) 

[XT I 



x J*(p, k)dpdk 



<ct\\w \\ B M\J\U + \\J\K b y 



(3.15) 



for all e E (0, 1], t > t , W G B a and J E A! A D B a , b . 

The proof of the above theorem is contained in Section 16.21 The 
remaining part of the argument follows the argument of Section 13.11 



4. Proofs of Theorems 13.21 and 13.51 
4.1. Proof of Theorem I3J2J. From f!2.45p and 02.441) we conclude that 



rl+7 



■(W e (t),J)-{U(t),J) 



< 



rl+7 



-W t (0)-U ,J(t) 



.1+7 



(Tl e {s),J(t-s))ds 



WeE (0,1], 



with TZ € (t) given by ( I2.37p . Estimates of the last term on the right 
hand side above shall be done separately for each term appearing on 
the right hand side of (12.371) . 
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4.1.1. Terms corresponding to TZ e ■ Denote 

f l+7_ f 1+ "/_ 

8 {t \t,k) := —W e (t,0,k) and G {t \t,k) := —Y e (t, 0, k). 
Lemma 4.1. For a given a G (0, 1] there exists C > such that 

/or aZZee (0,1], t > 0. 

Proof. Denote the expression on the left hand side of (14. ip by S a (t). 
From (EMD and fEOTp we conclude that 

/dk 
^ [\C(S^(t))(k)\ + |£(HegW(*))(fc)|] (4.2) 

Since \k\~ 2a R(k, k') is bounded when a G (0, 1] we can bound the right 
hand side of flt^D , with the help of f[2T2T]) , by C/C 7 and flO) follows. □ 
Let Auj e (p,k) := pu'(k) — 5 e u(p,k). For a given q G R define by 
mod(g, 1/2) the unique r G [—1/2,1/2) such that q = £ + r, where 
£ G Z. Divide the cylinder R x T into two domains: C - described 
below - and its complement C c . The first domain consists of those 
(p,k) for which either |mod(A; ± ep/2, 1/2)| > 1/4, or both points 
mod (A; ± ep/2, 1/2) belong to an interval [—1/2,0], or they belong to 
[0, 1/2]. We can write then 

e l+7 

— (AlJ e We(s), J{t - S)) =h + J 2 , 

where the terms on the right hand side correspond to the integration 
over the aforementioned domains. One can easily verify that there 
exists C > such that \Aou e (p,k)\ < Ce\p\ for (p,k) G C. We can 
estimate therefore 

\h\<Ce\\J\\ A )C T (4.3) 

On the other hand, when (p, k) G C c , with no loss of generality assume 
that 1/2 > k + ep/2 > > k - ep/2 > -1/2 (the other cases can 
be handled analogously) p is positive and k G (—ep/2, ep/2). Since 
\Au} e (p, k)\ < C\p\, for a given a G (0, 1] we can write 

\Aco t (p,k)\ Y\ | sin 

£T=±1 



d£g(t) 

dt 



7T 



aep\- 
2 J. 



< Ce 2a \p\ 2a+l 
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for some constant C > 0. Using the above estimate we obtain 
\h\<Ce 2a [ \p\ 2a+1 E\J(p,K(t-s,k))\ 

Jrxt 



r l+ 7 



n 



sin 



=±i 



7T ( k 



crep\ 



By virtue of (14.11) we get that 

|J 2 | < Ce 2o (s/C T + /C a;7 ;||j| U / a+ 

Taking the Taylor expansion of -R e (p, fc, fc'), up to terms of order e, it 
is also straightforward to conclude that 



(4.4) 
dpdk. 

(4.5) 



rl+7 



((C^-C)W t (s),J(t-s)) 



< e/C 



Summarizing, (14. 3p . (14. 5 p and (I4.6P together imply 



[ (TZ^\s),J(t-s))ds 
Jo 



(4.6) 



(4.7) 



< CeplU^t + Ce'H^lC, + /C a/y )|| J\\ , 



(2) 

4.1.2. Terms corresponding to1Z e . Straightforward computations, tak- 
ing the Taylor expansions of 0(k — ep/2) and R e (p, k, k') up to e, show 
that 

,1+7 pt pl+7 pi 

— J (K?\ S ),J(t-s))ds-—J (C(ReY e (s)),J(t-s))ds 

< CepW^t. 
From (I2.39P we conclude the following estimate. 
Lemma 4.2. Suppose that <p t : R x T — > R is such that 



1 * := sup 

e,p,k U e {p,k) 

Then, there exists C > such that 

rl+7 ^ 



< +oo. 



<CeK,,{\\J\\ A ,t 



(y e ( s )0 e ,j(t- s ))d s 

/or all J e S, e e (0,1], t > 0. 

Proof. The left hand side of (14. 9p can be rewritten as 



.4' 



(4- 



(4.9) 



rl+7 /•* 



(w e F e (s)r e , J(t-s))ds 



(4.10) 
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where F e (p,k) := (p e (p,k)u e 1 (p,k). Using ( 12. 39ft we can estimate this 
expression by 



.1+7 /■* 



/ (d s Y e (s)T e , J(t-s))ds 
Jo 





e i+ 7 /• 


+ e 






2 7o 



(w e (s)r £ , J(t - s))ds 

(4.11 

Thanks to (12.211) and (14. 8p the second term is bounded by el?* || JW^JC^t. 
On the other hand, integration by parts allows us to estimate the first 
one by 



r l+7 ft 



(Y e {s)T e ,{iB + C)J{t - s))ds 



(4.12) 



£ l+7 



(Y e (s)T e ,J(t-s)) 



s=0 



<C/C 7 e(*||J|U i + ||J|U), 



by virtue of (1231]) and (Oil . 
Since 

sup — — - — — — < +oo 

e,fc,fc',p 0J e (p,k) 

from the above lemma we conclude that 



□ 



(4.13) 



,1+7 r l 



(C(ReY e (s)),J{t- s))ds 



<C/C 7 e(t||J|U ; + ||J|UO- (4-14) 



This ends the proof of (13.71) . 



□ 



4.2. Proof of Theorem 13.51 We maintain the notation from the ar- 
gument made in the previous section. Estimate (I4.4p can be improved 
since in this case there exists C > such that \Au e (p, k)\ < Ce\p\ for 
all (p, k), e G (0, 1]. Thus, we can write 

\h\ <C7C 7 e/ \p\E\J(p,K(t-s,k))\dpdk < C/C 7 e||J|U ; (4.15) 

Jrxt 

for all t > 0. With this improvement in mind we conclude that there 
is a constant C > such that 



r l+ 7 f t 



(K£\s), J(t-s))di 



< C€\\ J^A'jCryt 



(4.16) 



for all t > 0, e G (0,1]. Repeating the argument from the proof of 
Theorem 13. 21 for the term corresponding to 1Z € we conclude (13.131) . □ 
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5. Convergence rate for a characteristic function of an 
additive functional of a markov process 

5.1. Markov chains. Suppose that n > 0} is a Markov chain 
taking values in a Polish metric space (E, d). Assume that n - the law 
of £o - is an invariant and ergodic probability measure for the chain. 
The transition operator satisfies: 

Condition 5.1. (spectral gap condition): 

a := sup {IIP/Hi^) : \\f\\v M = 1, / JL l} < 1. 

Since P is also a contraction in L 1 ^) and L°°(7i) we conclude, via 
Riesz-Thorin interpolation theorem, that for any f3 G [1, +oo): 



11*711 



Z/(tt) 



for all / G Lq(tt) - the subspace of Lq{tt) consisting of functions satis- 
fying / fdn = 0, with k{P) := 1 — 1 2//5 — 1| > 0. Thus, Q N := En= pn 
satisfies 

WQNfWvw < (1 - a^y'Wfh^, V/ G LJ(tt). (5.2) 

Furthermore, we assume the following regularity property of transition 
of probabilities: 

Condition 5.2. (existence of bounded probability densities w.r.t. it) 
transition probability is of the form P(w, dv) = p(w, v)ir(dv), where the 
kernel p(-, ■) belongs to L°°{tt £g> 7r). 

5.2. Convergence of additive functionals. Suppose that ^ : E — > 
R satisfies the tail estimate 



1/(71-)) 



(5.1) 



7r(|$|>A)<-^ 



for some C > and all A > 1 and 



0. 



We wish to describe the behavior of tail probabilities P 
when k > for the scaled partial sum process 

! [Nt] 

z^ = jv^E*(W, *>o. 



(5.3) 

(5.4) 

> N K 

(5.5) 



n=0 



To that end we represent z\ N ^ as a sum of an L' 3 integrable martin- 
gale for f3 G [1, a) and a boundary term vanishing with N. Let \ be 
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the unique solution, belonging to Lq(it) for f3 G [l,a), of the Poisson 
equation 

X ~ PX = (5.6) 

In fact using Condition 15.21 we conclude that P\ G L°°(tt). There- 
fore the tails of x an d ^ under w are identical. We introduce an L 13 
integrable martingale letting: Mq := 0, 

N 

M N := Z n, where Z n := X (£n) - P X (Cn-i), N > 1 (5.7) 



n=l 



and the respective partial sum process m[ n ^ := iV _1 / a Mnvt], t > 0. 

Using the dual version of Burkholder inequality for L 13 integrable 
martingales, when (3 G (1,2), see Corollary 4. 22, p. 101 of [2D] (and 
also HU) we conclude that there exists C > such that 



\M N fy /IB <CN 1/I3 , VN>1. 



(5.1 



Lemma 5.3. Under the assumptions (15. 3p and f l5.4[) /or any k > 

and 5 G (0, «k) t/iere exzsi C > swca £/ia£ 



3 



(TV) 



> TV* 



Proof. Choose /3 G [1,«). We can write 

[iVt] 



(5.9) 



(5.10) 



n=0 



From (15.101) and Chebyshev's inequality we can estimate the left hand 
side of fl52D by 



Mi 



[JVi] 



> 



Mf[Wl > 



3 



< P 



x(£o)| > 

Mr 



l[Nt]\ 



> 



3 



+ 



c 



]yP(l/a+K) 



for some C > 0. On the other hand 



M, 



[Nt] 



> 



< 



E Mi 



[Aft] 



< 



Ct 



JS[P{l/a+K)- 



The last inequality follows from (15. 8p . Choosing /3 sufficiently close to 
a we conclude the assertion of the lemma. □ 
Suppose furthermore that an observable $ : E — > R is such that 
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Condition 5.4. J ^dir = and there exist a £ (1, 2), a.\ £ (0, 2 — a) 

and nonnegative constants c£,c~ , C* > such that c+ + c~ > and 



v 7 A Q 



tt(^ < -A) - ^ 
v 7 A a 



< 



(5.11) 



/or all X > 1. 

Let £ > 0} be an a-stable process with the Levy exponent 



iXp 



ip(p) := a 

where 

c*(A) := 

In what follows we shall prove the following. 



c # , when A < 0, 
cf, when A > 0. 



(5.12) 



(5.13) 



Theorem 5.5. Under the assumptions made above, for any 5 £ (0, a/ (a+ 

1)) t/jere exzsi C > snc/i £/ja£ 



Ee 



^ t (JV) _ „-*0(p) 



<c(i + w) 5 (t+i; 



w i (5 ' 14) 



/or allp £ R, i > 0, N > 1. 



Proof. In the first step we replace an additive functional of the chain 
by a martingale partial sum process. Using (I5.10P we conclude that 
there exists C > such that 



Ee 



ipZ 



(N) 



Ee ipM t 



(N) 



< 



c\ P \ 

N't/a 



Vi > 0, N> l,p £ R 



and (I5.14p shall follow as soon as we can show that for any 5 £ 
(0, a/(a + 1)) there exist C > such that 



Ee i P Af t (JV) _ e -#(p) 



<C(1 + H) 5 (t + 1) 



1 



(5.15) 



for all p £ R, £ > 0, N > 1. The remaining part of the argument is 
therefore devoted to the proof of (I5.15p . Denote Z^ >n := N~ 1 ' a Z n with 
Z n defined in (15. 7p and x the solution of Poisson equation (I5.6P with 
the right hand side equal to Introduce also 



h p (x) 



e ipx — 1 — ipx, 
Eh p (Zn,i) 



(5.16) 



and V (JV) (p) := -Nh p N] . Observe that Re /if < 0. 
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Lemma 5.6. There exists a constant C > such that 



|^W(p)-^(p)|< " |p|(l + bl) 



C 



(5.17) 



In addition, for any bounded set A C R and (3 G [l,a) there exists 
C > such that 



N I sup 

AeA 



h, 



V(v) + X 



ir(dv) < C\pf(l + \p\), VpGR, N>1. 

(5.18) 



Proof. Denote ^f(w,v) := ^(u) + Pxi v ) ~ Px( w )- With this notation 
the expression on the left hand side of (I5.17P can be rewritten as 



N 



ExE 



h p [ N i/ a I P(w,v)Tr(dw)n(dv) -a h^—^dX 



c*(A) 



< 



N 



h, 



ExE 

N [ h p (A) 



N l l a 



p(w, v)ir(dw)7r(dv) 

(5.19) 



Here -F/v(A) := 7r( 1 J/ < iV^A). The first term on the right hand side 
can be estimated by 



CN i-ih 



sup 

E AeA 



hi 



m(v) + A 



ir(dv) 



where A is a bounded interval containing all possible values of Px{ v ) — 
Px(w). This expression can be further estimated by 

Cp 2 N l-2/a f (^(^l + l)^). 



As for the second term on the right hand side of (I5.19p . using integra- 
tion by parts, we conclude that it can be bounded by 



£T=± 



+00 



\KW\ 



< C*\p\N~ ai/a 



7r(a^ > iV 1/o A) 

oo lipA _ -y 



dX 



NX a 

dX < C\p\(l+ \p\)N^' a . 



(5.20) 



The first estimate follows from ( 15. lip , while the last one follows upon 
the change of variables A' := Xp and the fact that 1 < a + a\ < 2. 
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Concerning (I5.18P expression appearing there can be estimated by 



sup 

r AeA 



\h p (\' + \N- 1/a ) 



N / \h p (\')fF N (d\') + f3N 1 ^ a 
Jr 

x \h' p (A' + XN~ 1/a )\F N {dX'). 
The first term is estimated by 

r+oo 

(3N / \hp (A'/" 1 \h' (A')k(|*| > N^XyX' < C\p\ 



(5.21) 



for some C > and all p G R, N > 1. Since A is bounded the second 
term on the other hand is smaller than 

CWi-V" A suv \h p (X' + XN^/^f" 1 \h' p (X' + XN^ 1/a )\F N (dX') 

Jr AgA 



(5.22) 



Hence, fl5TT8l) follows. □ 
In fact, in light of the above lemma to prove the theorem it suffices 
only to show that for any 5 G (0, a/ (a + 1)) one choose C > so that 



Ee 



(JV) 



<^(i + bi) 5 (t + i; 



(5.23) 



for all p G R, t > 0, N > 1. To shorten the notation we let M^jv := 
Mj/N 1/a . Since {M„, n > 0} is adapted and E[Z N>n+1 \ F n \ = 0, we 
can write 



e ipZ Ntj+1 | jr. 



e^-Wl + E 



hp (Zjvj+i) | J 7 
To derive a recursive formula for 

Wj := exp{^ (JV) (p)i/^}Eexp{ipM iiJV }, 

write 

W j+1 - = eM^ (N) (p)(j + l)/N}E{^ M ^[h p (Z N;j+1 ) - h^} 
+ exp{^ N \p)(j + l)/iV}(l + h p N) - e h ? N) )Ee ipM ^ N . 
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Since Mq = 0, adding up from j = up to [Nt] — 1 and then dividing 
both sides of obtained equality by exp-fyW (p)[Nt]/N} we obtain that 

E[exp{zpM t (A °}] -exp{-tP iN) (p)[Nt}/N} (5.24) 

[Nt]-1 [Nt]-1 
3=0 j=0 



.r 



where 



e NJ := exp{^ (A °(p)(j + 1 - [Nt])/N}e ipM ^ N . 



We denote the terms appearing on the right hand side of (15.241) by / 
and II and examine each of them separately. As far as 77 is concerned 
we bound its absolute value by 



Nt 



i + m - e h - 



(N) 



(5.25) 



and since Re hp < we obtain the following 



\II\ < Nt 



< 



Ct\p\\l + |p|) ; 
N 



(5.26) 



for some C > 0. 

Fix K > 1, to be adjusted later on, and divide the set = 
{0,...,[Nt} - 1} in t = [[Nt]/K] + 1 contiguous subintervals, t of 
size K and the last one of size K' < K, i.e. 



A N = [J X m , I m Dl n = for m ^ n , 

m=l 

1-m = {jm, ...,j m + K -1}, m = !,...,£-! 



and 



1-1 = {jm, ■ ■ ■ , jm + K'} 



with K' < K. Here [a] stands for the integer part of a G R. To simplify 
the notation we shall assume that K' = K. This assumption does not 
influence the asymptotics. 

We need to estimate the absolute value of 



1 = E J2^{eNAh P (Z N , J+ i) ~ h p N) ]} =h + h 



(5.27) 
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k=\ jei k 



h:= t E { e ^EW%+i)-^] 
The conditional expectation in the formula for I\ equals 



h, 



7r(dv). 



The supremum of its absolute value can be estimated by 



2|b(-r)ll 



oo / SUp 

AeA 



h, 



ir(dv) < 



C\p\(1 + \p\) 
N 



for some constant C > 0. Here A is a bounded set containing and 
all possible values of Px{ w ) ~ Px( z ) f° r z,w E E. The last inequality 
follows from Lemma 15.61 On the other hand, since the real part of 
logejvj is non-positive we have 



Z\e N>j - e N , Jk \ < \p\ [KElh^l + E 



sup \M hN - M jkiN \ 



V j e l k 



According to (15.181) the first term in parentheses can be estimated by 
C\p\(l + \p\)K/N. Choose (3 G (l,a). From (EUD and Doob's inequality 
we can estimate the second term by CK 1 ^ /N 1 ' 01 . Summarizing we 
have shown that 



|/i| < C(t+1) (1 + |p|) 



, K 1 ^ K 



N x l c 



On the other hand, 

i 

I 2 = E E { e N,j k E E [ h p( Z N,j 



+1, 



N 



(5.2* 



T 



J l: 



k=i jei k 
e k-i 

k=l j=0 



where 



9n(w) :-- 



X(v) - Px{w) 
N l l a 



p(w,v)ir(dv) — 



\9n\\l^(w) 
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Fix e (1, a). Since / g N dn = and £ = [[Nt]/K] + 1, from Q52) we 
conclude that 

| J 2 | < (1 - fl-W)" 1 (^ + 1 

ESI , , /v 1 ^ 1 // 3 

< C(l - a^))- 1 ^ + l)|p|(l + |p|)V/»__. 

Hence, we have shown that 

|/| < C(t + !)(!+ | P ir + ^ + — 

Choose K = N 5 with 5 e (0, 1). It is easy to see that 

l-5>--l V/3e(l,a) 
a p 

thus the middle term on the right hand side is smaller than the first 
one. The optimal rate is obtained therefore when 

1 S r 1 1 

-- 15 = S + -o -1, 

a p p 

or equivalently 

l/a + 2 

IZsTT 

From the above we get, that the optimal rate of convergence is obtained 
when (3 is as close to a as possible, and for each Si < a /{a + 1) we can 
choose then C > so that 

Taking into account the above and (15.261) we conclude (I5.23p . □ 

5.3. Convergence rates in the central limit theorem regime. 

We maintain assumptions made about the chain in the previous section. 
This time however we assume instead of Condition (15.41) the following. 

Condition 5.7. / ^dn = and * G L 2 (n). 

We define the martingale approximation, using (15. 7p from the pre- 
vious section. Decomposition (I5.10P remains in force. Instead of ( 15.81) 
we have then an inequality 

(E\M N \ 2 ) 1/2 <CN^ 2 , ViV>l, (5.29) 
where C > is independent of N > 1. Let 

a 2 := EM 2 and rj)(p) := a 2 p 2 /2. 
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Define the partial sum process z\ by (15.51) with a replaced by 2. 
Our main result concerning the convergence of characteristic functions 
reads as follows. 



Theorem 5.8. There exist C > such that 



Ee 



ipZ, 



(N) 



-ti>(p) 



<C(l + b|) 



L t + 1 

7VV3 



(5.30) 



for dip E R, t > 0, N > 1. 

Proof. Denote Z N ^ n := N~ l / 2 Z n with Z n defined in (15. 1\ and Mj jN := 
Mj/N 1 / 2 . Using this notation as well as the notation from the previous 
section we can write that 



Ee 



ipZ 1 



(N) 



Ee 



ipM 



(JV) 



< 



c\ P \ 



Vt > 0, N > l, p e R. 



It suffices therefore to prove 



Ee 



ipMl 



(N) 



-t*l>(p) 



<c(i + b|> 



[ t + 1 

iVV3 



(5.31; 



Let /i p (x) be given by ( 15 .160 . We denote /i£ := — (ap) 2 /(2N) and 
W,- := exp{(ap) 2 j7(2iV)}E[exp{^M^}], 
Analogously to what has been done in the previous section we obtain 

W j+1 - W j = exp{(a P ) 2 ((j + l)/(2N))}E{e i ^[h p (Z Nij+1 ) - h^} 
+ exp{(ap) 2 ((j + l)/(2JV))}fl + - e^W^-" . 



(5.32) 



Adding up from j = up to [iV£] — 1 and then dividing both sides of 
obtained equality by exp{(ap) 2 [Nt]/ (2N)} we obtain that 

E[exp{ipM^ N) }] - exp{-(ap) 2 [Nt}/(2N)} 

[Nt] [Nt] 



3=0 



3=0 



where 



e NJ := exp{(<rp) 2 (j + 1 - [iVt])/(2iV)}e^/ Ar 



h (N) 
e p 



1/2 



Eejvj, 



Denote the absolute value of the terms appearing on the right hand 
side by I and 77 respectively. We can easily estimate 

for some constant C > and alH > 0, N > 1 and p6R. 

To estimate 7 we invoke the block argument from the previous sec- 
tion. Fix K > 1 and divide the set = {0, . . . , TV - 1} in I = 
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[[iVi]/.ftT] + 1 contiguous subintervals, I of size K and the last one of 
size K' < K. In fact for simplicity sake we just assume that all inter- 
vals have the same size and maintain the notation introduced in the 
previous section. Estimate (I5.27P remains in force with the obvious ad- 
justments needed for a = 2. Repeating the argument leading to (I5.28P 
with ( I5.29P used in place of (15. 8p we conclude that 

M<&(t + 1) + (5.33) 



On the other hand 

|/ 2 | < (i - a y l ( ^ + i ) < c(i - a 



"0H 



2 



K 

Choosing K = N 1 ^ 3 in the above estimates we conclude (I5.3ip . 

5.4. Convergence rates for additive functionals of jump pro- 
cesses. Assume that {r n , n > 0} are i.i.d. exponentially distributed 
random variables with Etq = 1 that are independent of the Markov 
chain {£ n , n > 0} considered in the previous section. Suppose further- 
more that V, 9 are Borel measurable functions on E such that: 

Condition 5.9. For some t* > 

6(w) >U, V weE, (5.34) 

there exist C* > and > 1 such that 

tt(0>A)<-^, VA>1. (5.35) 

We shall also assume that either ^/(w) := V{w)6{w) satisfies Condi- 
tion I5.4[ or it satisfies Condition 15.71 

Let t n := Xlfc=o 8(£k)i~k and let {X t , t > 0} be a jump process given 
by X t := £ nt where n t := n for t n < t < t n+1 . The process conditioned 
on £o = w shall be denoted by X t (w). The corresponding chain shall 
be called a skeleton, while 8~ 1 (w) is the jump rate at w. Note that 
measure 7t(dw) := 0~ 1 6(w)ir(dw) is invariant under the process. Here 

9 := E0(f o ). 

We also define processes 

Y {N) ._ N -x/p / V (x 8 )ds, (5.36) 

with j3 — a when this condition holds and j3 = 2, in case Condition 15.71 
is in place. 
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Observe that the integral above can be written as a random length 
sum formed over a Markov chain. More precisely, when (15. lip holds 

y t (N) = J^- a E *(4) + *f\ (5-37) 

where 

^ (7V); = jy^CU.X* (5.38) 

{£n, w > 0} is an E x (0, +oo)-valued Markov chain given by £ n : = 
(£n, T n ) and \P(u7,r) := fy(w)r. Its invariant measure n is given by 
fr(dw,dT) := e~ T 7i(dw)dr and the transition probability kernel equals 
P(w,T,dv,dr') := p(w,v)e~ T ir(dv)dT' . It is elementary to verify that 
this chain and the observable ^ satisfy Conditions I5.lti5.4l with the 
constants appearing in Condition 15.41 given by cf := T(a + l)cf, where 
T(-) is the Euler gamma function. In case Condition 15.71 holds we can 
still write floTBTI) and fl5~38D with a = 2. 

Consider the stable process {Z t , t > 0} whose Levy exponent is given 
by (15.121) with c*(A) replaced by 

7, m ._ / 0T a T(a + l)c~, when A < 0, . , 

c*{A) .- | Q-a V ^ a + 1)c + ; when \ > q l^ J J 

Let {B t , t > 0} be a zero mean Brownian motion whose variance equals 

c 2 := 20"V. (5.40) 
The aim of this section is to prove the following. 

Theorem 5.10. In addition to the assumptions made in Section [3771 
suppose that Condition \5.9i holds. Then, for any 

< 5 < <5* := min[a/(a + 1), (a 2 - l)/(aa 2 + 1)] (5.41) 

i/iere exists C > suc/i t/iai 

1 

(5.42) 

/or aWp e R, i > 0, JV > 1. 

7/j on i/ie oi/ier /jand t/je assumptions made in Section [TOI and Con- 
dition \5.9\ hold then for any 

0<5<5* := (a 2 -l)/(l + 2a 2 ) (5.43) 
tnere exists C > swc/i tnat 

Eexp{z P y t (iV) } - Eexp{zpi? t }| < C(|p| + l) 4 (t + 1) (-L + A 

(5.44) 



Eexp|^y t (JV) J -Eexp{zpZj < C(|p| + l) 5 (t + l) ^ 
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for allp G R, t > 0, N > 1. 

The proof of this result is carried out below. It relies on Theorem 
15.51 The principal difficulty is that definition of y} N ^ involves a random 
sum instead of a deterministic one as in the previous section. This shall 
be resolved by replacing nm appearing in the upper limit of the sum 
in (I5.37P by the deterministic limit [Nt], where N := N/6, that can 
be done with a large probability, according to the lemma formulated 
below. 

5.4.1. From a random to deterministic sum. 

Lemma 5.11. Suppose that k > 0. Then, for any 5 G (0, ci^k) there 
exists C > such that 



P [\n Nt - [Nt}\ > iV K+1 /« 2 ] < C ^ s 1 \ ViV > 1, t > 0. 
Proof. Denote 



(5.45) 



A + N := [n [Nt] - [Nt] > N 1 '^] , (5.46) 
A~ N := [n m - [Nt] < -N^ +K ] , A N := A + N U A~ N . 
Let k! G (0, k) be arbitrary and 

"[7Vt]+Af 1 /"2+K-l 

C N := m> Nl/a2+K ' 

n=[Nt] 

We adopt the convention of summing up to the largest integer smaller 
than, or equal to the upper limit of summation. Note that on A^ we 
have 

[JVt]+A rl / Q 2+«_i 

[Nt] — t[ftq > t n[Nt] — t[ftq > t* r n . 

n=[Nt] 

Furthermore 



r N :-- 



[Nt] 



[Nt] 



Hence, 



P [A+ n C N ] < P 1 1 [Nt] - t m I > t*N l ^ +K ' 
1 



< P 



[Nt] 

J2m n )rn-e] 

n=0 



> t*N K - 6r N 
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To estimate the probability appearing on the utmost right hand side 
we apply Lemma T5. 31 for the Markov chain {£ n , n > 0} and ^(w, r) : = 
6(w)t — 9. We conclude therefore that for any S e (0, c^ft') there exists 
C > such that 



p [A+ n c N ] 



< 



C{t 



N 5 



Vt > 0, N > I. 



(5.47) 



Since Etq = 1 and k' G (0, k) for any x G (0, 1) we can find C > such 
that 



P [C e N ] < P 



JV 1 / Q! 2+ K — 1 



< CP 



N 1 / 



Ari/a 2 

ATI - . 



r n < 



1 



n=0 

JV 1 /«2+ K - 1 



(5.48) 



< x 



n=0 



< C exp { -iV 1/a2+K /(x) } , V N > 1 



where J(x) := —(1 — x — lnx). The last inequality follows from the 
large deviations estimate of Cramer, see e.g. Theorem 2.2.3 of [TO] . 
Using this and (I5.47P we get 

p[i+]<p[ijnc,]+p[^]<p. 

Probability P [A^] can be estimated in similar way. Instead of Cn we 
consider the event 



C 



N 



[Nt] 

„=[7Vt]-Ar 1 /«2+ K + l 



and carry out similar estimates to the ones done before. 
5.4.2. Proof of (15 .42 \\ . Choose any k > 0. We can write 



E exp 



Eexp {ipZ t } 



□ 



(5.49) 



< 


E 


exp < 




> — exp < 




\,A N 


+ 


E 


exp < 


\ipYr) 


> — exp < 







E exp ^ipZ[ N ^ | — E exp {ipZ t } 

with An defined in (15.461) . The last term on the right hand side can be 
estimated by the expression appearing on the right hand side of (I5.42p . 
by virtue of Theorem 15.51 
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The first term on the right hand side can be estimated by C(t+1)N s 
for some 5 G (0, a 2 n) and C > 0. The second term is less than, or equal 
to 

\p\(El N + EJ N ), 

where 



1 A' 



N~ 1/a max 

me[[Nt]-N 1 / a 2+ K ,[Nt]+N 1 / a 2+'<-]nz 



[Nt] 
k=m 



and 



J N ■= N~ 1/a max 

me[[Nt]-N 1 / a 2+ K ,[Nt]+N 1 / a 2+ K ]nz 



(5.50) 



(5.51) 



Lemma 5.12. Suppose that k G (0, 1 — l/a^). Then, for any 5 G 
(0, a _1 (l — k — a^ 1 )) there exists C > such that 



El N < 



and 



C 



c_ 



EJ;v<^, ViV>l. 



(5.52) 



(5.53) 



Proof. First we prove (15.531) . Since ^(£o) is L 13 integrable we can write 
for any /3 G (1, a) 



EJ N < 



< 



E max 

m£[[Nt] -N 1 / a 2+ K ,[Nt]+N 1 / a 2+ K ]nN 



VP 



V/3 



(5.54) 



E 

m£[[Nt]-N 1 / a 2+ K ,[Nt]+N 1 / a 2+ K ]r\N 



Choosing (3 sufficiently close to a we conclude (I5.53p . 

Now we prove (]5.52p . Again we can use martingale decomposition 



m— 1 



m—l. 



Mr, 



n=0 



where 



m— 1 



n=l 



and x(') is unique, zero mean, solution of x ~ Px = ^> with P the 
transition operator for the chain {£ n , n > 0}. Using stationarity we 
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can bound 



El N < 



max 



N 1/a me[0,2JVV«2+K] nz 



m— 1 



2 



N l / a 



n=0 



max 



< — —j- E max |M„ 

N 1 '" mG[0,2Af 1 /"2+«]nz 



me[0,2Ar 1 /a2+«]nz 



Denote the terms on the right hand side by / 



(0 

N > 



1,2,3 respectively. 



One can easily estimate iffi < CN 1 / 0i . Also, to bound Ijp we can 
repeat estimates made in (I5.54p . as x is a l so ^ integrable and obtain 



r(3) 



for some C > 0. Finally, to deal with observe that by Doob's 
inequality for (3 £ (1, en) 



max 

^G[0,2^^ 1 /a2+' t ]nz 



IMJ < 



c 



E |M, 



2JV 1 /«2- 



|/3| 1//? 



We use again (15 .8p and conclude that 



E \M. 



2N 1 /°"2+ K I 



< CN {1/a2+K)/l3 . 



Summarizing from the above estimates we get 

/« < CN- S , 
where 5 is as in the statement of the lemma. 



□ 

Gathering the above results we have shown that for any k £ (0, 1 — 



ei- 



and 5i £ (0, a 2 K), S 2 £ (0, a 1 (1 — k — a 2 )) we have 



E exp 



[i P Y t 



(N)\ _ 



E exp < ipZ 



r(N) 



< C 



t 



Choosing k sufficiently close to (1 — ct 2 l )/ (aa 2 + 1) we obtain that for 
any S £ (0, (a 2 — l)(aa 2 + 1) _1 ) we can find a constant C > so that 



E exp < zpY^ 



.(TV) 



E exp < ij>Z 



AT)! 



(5.55) 



Thus, we conclude the proof of f)5.42p . 



5.4.3. Proof of (15.4=41) . In this case we can still write inequality ( 15.49D . 
With the help of Lemma 15.111 for any k > and 5 £ (0, a 2 n) we can 
find C > such that 

C(t+1] 



E exp 



|zpF t (7V) } - Eexp 



< 



N 5 



\p\(El N + Ej N ), 
(5.56) 
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where In, Jn are defined by (15.501) and (15.511) respectively, with a = 2. 
repeating the estimates made in the previous section we obtain that 

El N + E J N < CN 1 M"- 1+C * 1) 
for some C > 0. Using the above estimates and (15.301) we conclude 

tfESl). 

6. Proofs of Theorems 13.31 and 13.61 

6.1. Proof of Theorem Let N := e" 3 ^ 2 and J G A be a real 
valued function. Define 

i-Nt 

W N (t,p,k) := E ' 



W (p, K m (k)) exp |-^iV- 2 / 3 J uj'(K s (k))ds, 



(6.1) 

where {K t (k), t > 0} is the Markov jump process starting at k, in- 
troduced in Section 12.61 It can be easily verified that the Lebesgue 
measure on the torus is invariant and reversible for the process and we 
denote the respective stationary process by {K t , t > 0}. Its generator 
£ is a symmetric operator on L 2 (t) given by 

Cf(k) := f R(k,k')[f(k')-f(k))dk' 



z 



Here 



We also let 



Note that 



and 



-«(*)/(*) + | E V/GL 2 (T). 
*e{-i,i} 



tl (k) := -sin 4 (7rA;), e_i(ife) := 2 sin 2 (27r/c). 
3 



r(A;) := e_i(fc) + ei(fc). 



ti(k)dk = / t-\{k)dk 
Jz 



R(k) = jv{k) = 2sin 2 (vrfc) [l + 2 cos 2 (ttA;)] . 



4 

The process K t (k) is a jump Markov process of the type considered 
in Section 12. 6[ The mean jump time and the transition probability 
operator of the skeleton Markov chain n > 0} are given by 6(k) = 
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R-^k) and 



Pf(k) := 6{k) j R(k, k')f{k')dk' 



te{-i,i} 



e- t (fc) 
t(*) 



/eC(T), 



respectively. Probability measure ir(dk) = (l/2)v(k)dk is reversible 
under the dynamics of the chain. It is clear that Condition 15.21 holds. 
It has been shown in Section 3 of [13] that Condition 15.11 is satisfied. 

Let {Qt, t > 0} be the semigroup corresponding to the generator 
C. It can easily be argued that Q t is a contraction on L p (t) for any 
p G [1, +oo]. We shall need the following estimate. 



Theorem 6.1. For a given a e (0, 1] there exists C > 

C 



WQtfW 



Li(T) 



< 



B„, Vt>0, 



(6.2) 



;i + f)° 

for all f £ B a such that j T fdk = 0. 

The proof of this result shall be presented in Section I6~3"l We proceed 



first with its application in the proof of Theorem 
functional appearing in ( 16. ip shall be denoted by 

»JVt 



The additive 



Y} N \k) , 



AT2/3 



co'(K s (k))ds, 



or by y} N ^ in case it corresponds to the stationary process K t . It is of 

(6.3) 



the type considered in Section l5~4l with a = 3/2 and 

V(k) := u'(k)6(k). 
Since the dispersion relation satisfies 



u>{k) = \k\ 



d"(0) 



we have 



and 



with 



cu'(k) = sgn k 
tt(^ > A) 



d"(0) 



0(k 2 



0{k< 



1/2 



for fc«l. 



(6.4) 



1/2 



for k < 1 



A 3 / 2 



^+ . = 2 -i/4 3 -5/2 7r l/ 3 r a /// 0ll 3/4 
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and some C* > 0. Condition 15.41 is therefore satisfied with a = 3/2 
and arbitrary a.\ < 1/2. Since ^>(k) is odd we have c~ = c+. On the 
other hand, jump mean time 9{k) satisfies (I5.35P with a 2 = 3/2. 

We can apply to Y^ N ^ the conclusion of part 1) of Theorem 15. 101 In 
this case c*(A) = c, where 



2 l/2 



"°° sin 2 x 



x 



5/2 



dx. 



Since the integral on the right hand side equals 4^/^/3 we obtain (I3.5p . 
Let Z t be the corresponding symmetric 3/2-stable process. From (13.91) 
we obtain W(t,p) = W(p)Ee~ ipZt . Therefore, for any J(-) E S 



\(W N (t),J)-(W(t),J)\ 



J*(p, k)E W(p, K Nt (k))e~ ipY t N){k) - W(p)e- ipZ 



(6.5) 
dpdk . 



Let 1 > (3 > 1/3. The left hand side of (16.51) is estimated by E\ + 
Ei + E3, where 



E 1 :-- 
E 2 :-- 
E 3 :-- 



RXT 



r(p,k)E\w(p,K Nt (k)) 



ipY 



(N) 



r{p, k)EW p (K Nt (k))e t(i—N~@) dpdk 



RXT 



r(p,k)w(p)E 



e - ipY t"Z N -h {k) - e ~ ipZt 



dpdk 



and W p (k) := W(p, k) — W(p). The first term can be estimated as 
follows 



Ei < ||u; , || 00 tiV 1 / 3 - /3 f \p\\J(p,k)\E\W(p,K m (k))\dpdk 

J RXT 



< CtN 1/3 -P\\J\\ A[ \\W 



RXT 
.4 



(6.6) 



To estimate the second term note that by the Markov property we 
can write 



W p {K Nt {k))e 



-ipY 



(iV) 



t{l-N~P) 



(k) 



(6.7) 



Q m - H W P (K 



Nt(l-N-P) 



Term E 2 can be therefore estimated by 



(k))e 



-ipY 



(N) 



t(l-N-P) 



(fc) 



\J\\a> sup \\Qm-PtW p \\ L i^). 
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Invoking Theorem 16.11 we obtain 

sup \\Q N ,- H W P \W) < j^zfi \\W\\ Ba , Vt > 1, N > 1. 



PGR 



As a result we conclude immediately that 

C 



W\\ Ba \\J\U> 



(6.8) 



AT«(i-/3) 

To deal with E$ note that by reversibility of the process K t it equals 



W(p)E 



ipY 



(N) 



-ipZt 



dp 



where J{p) := f J(p, k)dk. We obtain that 



E 3 < 
+ 



W(p)EU*(p,K t(1 



(l-N-P)j 



-ipY 



(N) 



W(p)E 



J*(p)W(p)E 



-ipY 



ipY { 

t(l-2N~P) g 



(AT) 

t(l-2iV-/3) 



(N) 



-ipZt 



t(i-N-P) _ 

dp 
dp 



-ipY 



(JV) 



dp 



E3I + #32 + E- 



33 ■ 



From this point on handle this term similarly to what has been done 
before and obtain that 



E S i < 



Ct 



and 



#32 < 



c 



\\J\\a\\ w \U 



w\ 



(6.9) 



(6.10) 



provided b > 1. 

Term .E33 can be handled with the help of Theorem 15.101 with a = 
«2 = 3/2 and an arbitrary a.\ < 1/2. Therefore, for any 5 < 2/13 we 
can find a constant C > such that 



E 33 <^(t+l)\\W\\ A \\J\\Ai- 



(6.1i; 



We have reached therefore the conclusion of Theorem 13.31 with the 
exponent 7' as indicated in ( 13. 2p . 

□ 
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6.2. Proof of Theorem 13.61 The respective additive functional in 
this case equals 

1 f m 

Y t N) ( k )--=J^ u>(K s (k))ds, 



where N := e 27 . The observable \& is given by (16.31) . From (16 Ah we 
get 



d"(0)fc 



«(0) + ^A; 2 (1 + 0(A; 2 )) 



-1/2 



for fc < 1. 



and the asymptotics of the tails of ty(k) is given by 

tt(|*|>A)<^ 

for some C* > and all A > 1. The observable belongs therefore to 
L 2 (ir) and, since it is odd, its mean is 0. Conditions 15.71 and 15.91 are 
therefore fulfilled. The latter with «2 = 3. We also have 

c:=9a 2 , (6.12) 

where 

a 2 = j[ X \k)-{P X )\k)]n{dk) (6.13) 

and x is the unique zero mean solution of equation 

X-P X = V. (6.14) 

The result is then the consequence of the argument made in Section 
and Theorems 15.101 and 16.11 



6.3. Proof of Theorem 16.11 Denoting f t := Q t f we can write, using 
Duhamel's formula 

ft = Stf + \ Yl I S t - s t-,(k)(c,,f s )ds, (6.15) 

'£{-1,1} ° 

where 

S t f(k) := e-*Mf(k). 
Let H := [A G C : Re A > 0], 

r+00 

f(\):=(\-C)- 1 f= / e~ Xs f s ds, AGH, 

Jo 

fo(A) := -L- (6.16) 
r + A 
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and f t (A) := (/(A), e t ), t = ±1. Formula ( 16.160 extends to the resolvent 
set of the generator C in L 2 (t), that contains in particular C \ [— M, 0], 
with M : = (4/3)11^1100 + 1. 
From (16.150 we obtain that 

'«4(f) + E (6 - 17) 

v 7 ie{-i,i} ' 



Vector f T (A) := [f-i(A), fi(A)] satisfies therefore 



f (A) = 2T 1 



4A 

y 



4A 

y 



(6.18) 



where 



21(A) = 
a(A) := 1 



a 4 (A) := - 

g T (A) := [s_ 1 (A),g 1 (A)] and 

f(k)c L (k)dk 



a(A) O-i(A) 

ai(A) a(A) 

e_i(fc)ei(fc) 
T A + x(k) 

c 2 {k)dk 



dk. 



A + t(fc) 
Let A (A) := det 21(A) and 



Observe that 
and 



" AX): ~- A A + «(*)■ 
o(0) = -a_i(0) = -ai(0) 



x A + r(/c) ' 

t = ±1, A e c \ [-M,o]. 

z u {k)dk 



A(A) = A[A6_!(A)6i(A) + 6_ 1 (A)a 1 (A) + a„ 1 (A)6 1 (A)]. 



(6.19) 



(6.20) 



From (16.191) we get that A(0) = 0. In addition, from (16.200 we can see 
that D(X) := A(A)A~ 1 is analytic in C \ [-M, 0]. In addition, 



lim _D(X) = 6_i(0)ai(0) + a_i(0)&i(0) > 0. 

A^0,A6H 

Hence, there exist g > and c* > such that 

|A(A)| > c*|A|, VAgh, \\\<g. 



(6.21) 
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It can be straightforwardly argued that g can be further adjusted 
in such a way that ( I6.2ip holds on the boundary C of the rectangle 
(_ M ,0) x {-q,q). 

From (I6.18P we obtain that 

lW = A" 1 (^) n t (^) , t = ±l, (6.22) 

where 

n_i(A) := a(A)fl_i(A) - a_i(A)fli(A) 

and 

ni (A) := -a 1 (A)0_ 1 (A) + a(A)fli(A). 

Using the fact that J T /(i/c = 0, after a straightforward calculation, we 
obtain that 

n t (A) = A MA)so(A) - & t (A)&(A)] , t = ±1, (6.23) 

where 

flo(A) := / f Q (X,k)dk. 



Using the well known formula, see e.g. Chapter VII. 3. 6 of 

where JC is a contour enclosing the L 2 spectrum of £, that as we recall 
is contained in [— M, 0] . It is easy to see that 

|f t (A) | < C||/|| B JAr\ A g H, |A| < g, l = 0, ±1 (6.24) 

for an appropriate g > 0. This is clear for i = 0. Form this and (16 .21ft . 
( 16.23ft we conclude that (16.24)) holds also for i — ±1. Therefore, we 
can use as /C the boundary C of the rectangle, mentioned after ( 16 .21ft . 
oriented counter-clockwise. From ( 16 . 1 Tft we get 

/<= 7 e -'X(t + 4V3)A(4A/3)- (6 ' 25) 



Note that 



|5 t /|Ui (T )= / e-^«|/(*)|d*<|ll/lk, 



where C := sup x>0 x a e _a; . Consider the term 7(t) corresponding to 
the integral on the right hand side of (16.251) . We can write I(t) = 
Y2i=i Ii(t), where Ii(t) correspond to the sides of the rectangle {— M} x 
(-q,q), [-M,0] x {-g}, [-M,0] x {q}, {0} x {-g,g) appropriately 
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oriented. The estimations of i — 1, 2, 3 are quite straightforward 
and lead to the bounds 



Li(T) 



< 



c 



t + 1 



L i (T ), i — 1, 2, 3. 



(6.26) 



To deal with 7 4 (t) observe that, thanks to (16.231) . it equals to Ia(t) + 
/ 42 (t), where 



E ^f e!%jrl) (y)^ 

<.e{-i,i} e 



J = 1,2 



and 



a t (w) go [}v ) 



9i,2{v) 



-b L (iv) fl t {iv) 



x + iv x + iv 

The asymptotics of J 4j i for £ 3> 1 is, up to a term of order ||/||li( T )/£, 
the same as 

4i(*)-^ E ^ [ e lvt (Fg iA )(4zu/3)du 

t6{-l,l} 

for some C°° function F{iv) supported in (— p, p) and equal to D _1 (i^) 
in (— p/2,p/2). Denoting 



f)(A) := 



go (A) 
t+A 



we can write 



4iWIUmt)<^ E \^- idk 
/^'K)(f(-7))-^)(f 

f)(4i(i/ + 7r/t)/3) 



4 £ 



te{-i,i} 



47T 



^2 / e_ t dfc 



e{-i,ir T 



20 
2q 

2g 



[fj (4ii//3) - J) (4z(i/ + ?r/t)/3)] i 



(Fa,) 



Aiv 
IT 



for sufficiently large t (the support of F is contained in (— g, g)). The 
first term on the utmost right hand side can be estimated by 



C 

T 



2q q 

|flo (4z(i/ + 7r/t)/3)|di/< - 

■2g t 



\f\dudk C 
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As for the second term, it can be estimated by 



C 



2p 



* te ^ 1} A7-2pk + 4zz//3||r + 4z(z/ + 7r/t)/3| l ' j 

+ C ^ [ f 2p c_ L dkdv f \f\dk 



t f^, x JiJ-2 P |t + 4zV/3| J T |r + 4ii//3||r + 4i(z/ + 7r/t)/3| 



The first term is less than 

C ^ f i-M f \f\dk f 2p dv 



1 J^ )X }Ji tl+b J? xa J-2 P W\ 1 - b W + */t\ 1 - a 



(6.28) 



for some b G (0, 1/2). Since for any a, b > there exists C > such 
that 

r2p dv C 

rrrn i?— < ^ — Vx>0 (6.29) 

^\v\ l ~ b \v + x^- 11 ~ x 1 -^ 

expression in (16.281) can be estimated by C\\f\\ Ba /t a+b . Finally the 
second term in (16.27P can be estimated by 

C f f 2p dv f \f\dk C 



t J r J_ 2p \v\ x - a / 2 \i/ + 'ir/t\ 1 - a / 2 J T r a ~ t a 
by virtue of (16. 29ft . Summarizing, we have shown that 

C 



'4,H^||Li(T) < f t + iy \\J 

for some C > and all t > 0. The estimates for HZ^WHl^t) are quite 
analogous. 
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